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In this paper we use a spin kinetic equation to study spin polarization dynamics in ID wires 
and 2D channels. This approach is valid in both diffusive and ballistic spin transport regimes and, 
therefore, more general than the usual spin drift-diffusion equations. In particular, we demonstrate 
that in infinite ID wires with Rashba spin-orbit interaction the exponential spin relaxation decay can 
be modulated by an oscillating function. In the case of spin relaxation in finite length ID wires, it 
is shown that an initially homogeneous spin polarization spontaneously transforms into a persistent 
spin helix. An interesting sound waves echo-like behavior of initially localized spin polarization 
packet is found in finite length wires. We show that a propagating spin polarization profile reflects 
from a system boundary and returns back to its initial position similarly to the reflectance of sound 
waves from an obstacle. Green's function of spin kinetic equation is found for both finite and 
infinite ID systems. Moreover, we demonstrate explicitly that the spin relaxation in 2D channels 
with Rashba and Dresselhaus spin-orbit interactions of equal strength occurs similarly to that in ID 
wires of finite length. Finally, a simple transformation mapping ID spin kinetic equation into the 
Klein-Gordon equation with an imaginary mass is found thus establishing an interesting connection 
between semiconductor spintronics and relativistic quantum mechanics. 
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I. INTRODUCTION 

Dynamics of electron spin polarization in semiconduc- 
tor structures^ has attracted a lot of attention recently 
in the context of spintronics^. Typically, the most im- 
portant mechanism of spin relaxation in semiconductors 
lacking inversion sym met ry is the D'yakonov-Perel' spin 
relaxation mechanisnPS, intimately related to the spin 
splitting of the electronic states. Numerous studies of 
D'yakonov-Perel' spin relaxation in the past have dom- 
inantly concentrated on spin relaxation in infinite two- 
dimensional (2D) systems^E^l while the electron spin 
relaxation in confined geometries is not yet well under- 
stood. This problem, however, is of crucial importance 
since actual spin-based electronic devices^ are normally 
of a finite size. 

There are only several examples in the literature where 
the effects of boundary conditions on D'yakonov-Perel' 
spin relaxation have been explored theoretically and/or 
experimentally. In particular, investigations of spin re- 
laxation in systems with boundaries include spin re- 
laxation in 2D channelfPHUl, 2D half-s paceP , 2D sys- 
tems with antidotiP^, large quantum dot d 26 * 27 1 , and one- 
dimensional (ID) finite-length wires^ and rings^. The 
common conclusion of these studies is that in the diffusive 
regime of spin transport the introduction of boundary 
conditions results in an increased spin lifetime. In par- 
ticular, in Ref. 28 the present authors studied spin relax- 
ation in finite length quantum wires using drift-diffusion 
equations approach. It was demonstrated that a persis- 
tent spin helijpQ^ spontaneously emerges in course of 
relaxation of homogeneous spin polarization. 

In this paper, we extend the results of our previous 



study 2 ^ applying spin kinetic equation approach to the 
problem of spin relaxation in ID wires and 2D channels 
with spin-orbit interaction. The formalism of spin kinetic 
equation describes both diffusive and ballistic regimes 
of spin transport at any value of spin rotation angle 
per mean free path. Therefore, the spin kinetic equa- 
tion is mo re gene ral than traditional spin drift-diffusion 
equations ^ 8 * 32 * 33 ] and can be strictl y justified by using 
quasiclassical Green's function s] 19 * 34 !. Basically, the spin 
kinetic equation describes the spin polarization dynam- 
ics on shorter time and space scales than those of spin 
drift-diffusion equations. Consequently, the application 
of spin kinetic equation to the problem of spin relaxation 
provides new insights into the spin polarization dynam- 
ics. In particular, in this paper we report the bound- 
ary spin echo effect in which a localized spin polariza- 
tion profile reflects from a sample boundary and returns 
with a decreased amplitude to its initial position simi- 
larly to the usual spin echo of sounds waves. This type 
of spin echo is essentially different from the spin echo 
in nuclear magnetic resonance. The later is related to 
nuclear spins dephasing in an applied magnetic fields 
Moreover, we find a simple transformation that maps 
the spin kinetic equation into the Klein-Gordon equation 
with an imaginary mass, which is a relativistic analog of 
the Schrddingcr's equation. It is worth mentioning that 
the Klein-Gordon equation with an imaginary mass de- 
scribes tachyonS^KM . hypothesized particles that travel 
faster than light. Consequently, we suggest that semi- 
conductor spintronics structures have a potential to be 
used as a laboratory test bed for relativistic quantum 
mechanics. 

Our paper is organized as follows. In Sec. [IT] we intro- 
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duce a spin kinetic equation and derive a set of equations 
for spin polarization components and boundary condi- 
tions. We also demonstrate that the spin kinetic equa- 
tion can be transformed into the Klein-Gordon equa- 
tion. Next, in Sec. we employ the spin kinetic 
equation to study spin relaxation in infinite and finite 
length wires with Rashba spin-orbit interaction. In par- 
ticular, in Sec. |IIID| we present the boundary spin echo 
effect. The spin relaxation in two-dimensional channels 
with Rashba and DresselhauiP^ spin-orbit interactions of 
equal strength is investigated in Sec. |IV| We demonstrate 
that such a problem, for channels in a specific direction, 
can be mapped into the problem of spin relaxation in 
finite length wires with Rashba interaction only. The 
results of our investigations are summarized in Sec. |Vj 



II. THEORETICAL FRAMEWORK 

A. Spin kinetic equation 

The main goal of this paper is to study the kinetics 
of spin relaxation in wires made of 2D quantum well 
or heterostructure with Rashba spin-orbit interaction^. 
Therefore, we first introduce the Hamiltonian for an elec- 
tron in 2D space in the presence of spin-orbit interac- 
tion and all important parameters that will be used later 
in the spin kinetic equation. The standard Hamiltonian 
with the RashbaSH term is given by 



H — Hn + Hr — 



2m 



a (a X p) • z, 



(1) 



where p = {jp x ,p y ) is the 2D electron momentum oper- 
ator, m is the effective electron's mass, a is the Pauli- 
matrix vector, a is the spin-orbit coupling constant and 
z is a unit vector perpendicular to the confinement plane. 

It is not difficult to demonstrate^ that in the case of 
Hamiltonian (JlJ the quantum mechanical evolution of a 
spin of an electron with a momentum p can be reduced 
to a spin rotation with the angular velocity Q = 2ap/h 
about the axis determined by the unit vector n = pxz/p. 
In this way, the spin-orbit coupling constant a enters into 
equations through the parameter r\ = 2am.hr 1 , which 
gives the spin precession angle per unit length. 

Besides this evolution, 2D electrons experience dif- 
ferent bulk scattering events such as, for example, due 
to phonons or impurities. These scatterings random- 
ize the electron trajectories. Correspondingly, the di- 
rection of spin rotation becomes fluctuating what causes 
average spin relaxation (dephasing). This i s th e famous 
D'yakonov-Perel' spin relaxation mechanism!^ The time 
scale of the bulk scattering events can then be character- 
ized by a single rate parameter, the momentum relax- 
ation time r. It is connected to the mean free path by 
I = vt, where v = p/m is the mean electron velocity. 
To take into account these scatterings we use a kinetic 
model of spin transport presented below. When a char- 
acteristic system size L>f and the time scale of interest 



is much longer than r , the spin kinetic model yields the 
spin drift-diffusion equations (such as, e.g., reported in 
Ref. |33}. 

In the semi-classic approximation the kinetic equation 
for electron spin polarization can be written as (see, e.g., 

Refs. 023 uni H3I3J} 



d_ 
dt 



£ 

m 



V S r 



rip x Sp 



St{S p }, 



(2) 



where S p (r,i) is the vector of spin polarization of elec- 
trons moving with momentum p, and 5i{S p } is the colli- 
sion integral describing electron scattering processes. In 
the r-approximation the collision integral is given by 



St{S p } 



:(Sp-(S p )), 



(3) 



where the angle brackets denote averaging over direction 
of electron momentum. The conditions of applicability 
ofEq. Q can be found in Ref. [03 The collision integral 
([3]) corresponds to the elastic scattering of electrons by 
strong scatterers with a characteristic time r between the 
collisions. Note that the collision integral ^ conserves 
the total spin polarization redistributing spin polariza- 
tion between electrons moving in different directions. For 
ID case the average spin polarization simplifies to the fol- 
lowing expression (S p ) = (S+ + S~)/2, where S+ and S" 
are the spin polarizations of electrons moving along the 
wire in the positive (with momentum p = mve x ), and 
negative (p = — mve x ) directions with the average veloc- 
ity v. Thus, the kinetic equation ^ for ID wire takes 
the form of the system of two vector equations 



" vS-)s+ = -ne y xS+-^(S+-S-), (4) 



dt dx 
d 



2t 
I 



dt ^]S- = ^xS---(S--S^ 



(5) 



This system of equations should be complimented by ini- 
tial conditions for spin densities S + and S~ 



S+(M = 0) = S+(a;), 
S-(M = 0) = So(z), 



(0) 
(7) 



and by boundary conditions for finite length wires. The 
boundary conditions conserving the spin polarization in 
elastic scatterings at the boundary T = [x = 0, x = L] 
have the form 



(S+-S" 



= 0. 



(8) 



Taking the sum and difference of Eqs. (|4j |5|) we easily 
obtain 



dS dA n 

— +v— +ne y x A = 0, 
at ox 

OA dS n A 



(9) 
(10) 
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where the following notations are used: S = S + + S _ and 
A = S + — S — . As we are mainly interested in finding 
the total spin polarization S, A can be eliminated from 
Eqs. (55) via the following transformation. First of 
all, we multiply Eq. ( 10 ) by e'/ r and rewrite it as 



^+^g + n e -e»>S = 0. (11) 

Then, Eq. ^ is multiplied by e*/ r and is differentiated 
with respect to time 



0_ 
dt 



d d(e t / T A) 
dx dt 



+ fte y x 



afe'/ T A) 



dt 



= 0. 



Finally, we substitute d (e*/ r A) /dt from Eq. (IT] ) into 
Eq. (12 1. The resulting equations for spin polarization 



can be presented as 

d 2 S IdS _ v2 d 2 S _ 2Qve dS 
dt 2 t dt dx 2 y dx 

+ n 2 {S- S y e y ) = 0, (13) 

where S y is the y-component of S. 

Next, we would like to reformulate the boundary con- 
dition given by Eq. (8) in terms of the function S only. 
It is easy to notice that Eq. (8 ) corresponds to A| r = 



Substituting this boundary value of A into Eq. ( 10 ) we 
obtain 



v— — h ile v x S 
dx 



0. 



(14) 



Let us introduce a complex polarization 



(18) 



It is straightforward to show that Eq. (13 1, and boundary 



conditions ( 14 1 can be rewritten in a more compact form 
using S: 



d 2 S 1 as 
W + r~di 



d_ 

dx 

dS _ 

v— lilS 

dx 



m) S = 0, 
= 0. 



(19) 
(20) 



Moreover, it follows from Eqs. (19) and (20 1 that they 



(12) have only one stationary solution 



S = Soe* 



(21) 



where So is an arbitrary com plex constant. This solu- 
tion is so-called spin helisPD, which is persistent in ID 
Rashba wires despite electron collisions. Taking the real 
and imaginary parts of Eq. (21 ) with 5o = Ae lS , we find 



the usual representation of spin helix 

S x = Acos(i]x + S), S z — j4sin(?7a; + S), 



(22) 



where A is the spin helix amplitude, and 6 is its phase. 
Since the solution (22) of Eqs. (19) and (20) is the only 



stationary one, any initial spin polarization distribution 
eventually transforms into the spin helisP^, in some cases, 
however, having zero amplitude A — 0. 



Previously, the same form of boundary condition was de- 
rived using the Green's function methocP^. 

Moreover, since Eq. ( 13 ) is the second order differen- 



tial equation with respect to time, one must specify both 
the spin polarization and its time derivative at the initial 
moment of time t = 



S(z,t = 0) = S (z), 



as 

~dl 



= S (x). (15) 



Note that if we know A at the initial moment of time 
t — then we can find the first time derivative of S at 
t = using Eq. In particular, if A(x,t = 0) = 

then it follows from Eq. Q that S (x) = 0. 

Considering y-components of Eqs. (13)-(14) we find 



that S y is not coupled to any other component of spin 



polarization. Specifically, Eqs. ( 13 (-( 14 ) for S y can be 
rewritten as 



d 2 S y ^ IdS, 



d 2 S 1 

r dt dx 2 



dSy 

dx 



= 0. 



(16) 
(17) 



Consequently, selecting S y (x,t = 0) = we can safely 
take out S y from our consideration. 



B. Klein-Gordon equation 

We can exclude rotations of spin polarization vector 
that are still present in Eq. ( 19 ) and Eq. ( 20 ) by intro- 
ducing a complex field u via 



u(x, t) 



--S(x,t). 



(23) 



It can be shown that Eq. (19) and Eq. (20) transform 
into 



dt 2 V 'dx 2 4r 2 " °' 
du 
dx 



(24) 
(25) 



The transformation given by Eq. (23 1 is useful, in partic- 



ular, since it allows us writing down the Green's function 
Sg(x, t) of Eq. ( fl9| through the known Green's function 
of Eq. (24) (we will take advantage of this fact below). 
We note that Eq. (24) coincides with the well-known 
(see, for example, Ref. S3]), but 
(see the sign of the last term) . 



Klein-Gordon equation 
with an imaginary " mass 
In relativistic quantum mechanics such an equation de- 
scribes hypothetical particles called tachyongSSHSSl g e _ 
cause of the transformation (23), the tachyon physics 
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is somewhat 'hidden' in the dynamics of spin polariza- 
tion and can be recovered by an inverse transforma- 
tion. Analogies between certain exotic particles and con- 
densed matter systems have attracted a strong attention. 
For example, electrons in graphene^ and graphite^ ex- 
hibit properties of Dirac fermions, signatures of mag- 
netic monopoles were spotted in spin ice materials-^, 
and some exotic quasi-particle excitations in a variety 
of condensed-matter systems show a similarity with Ma- 
jorana fermiona^. Therefore, we believe that our ob- 
servation that is exciting in itself will stimulate further 
theoretical and experimental work in this area. 



III. SPIN RELAXATION IN WIRES 

In this section we present our studies of kinetics of 
spin relaxation in ID infinite and finite length wires with 
Rashba spin-orbit coupling using the spin kinetic equa- 
tion approach. Our analysis is essentially based on Eq. 



( 19 1 supplemented by initial and, where appropriate, by 
boundary conditions ( 20 ) . It is assumed that a thin 



narrow wire is made of a semiconductor heterostructure 
with spin-orbit coupling and the conduction electrons in 
the wire occupy the lowest size-quantization level cor- 
responding to transverse confinement. In addition, the 
length of finite length wires is considered to be much 
longer than the phase coherence length. Therefore, the 
electron transport in the direction along the wire is de- 
scribed in terms of the classical transport regime. 



Relaxation of homogeneous polarization in 
infinite wires 



Let us start the analysis of solutions of Eq. ( 19 ) with 



the most simple case, namely, with a problem of relax- 
ation of homogeneous initial spin polarization in an in- 
finite wire, — oo < x < oo. In this case, assuming that 
solutions of Eq. 



this equation anc 



19 1 do not depend on x, we can rewrite 



initial conditions as 



d 2 ^ 
dt 2 



IdS 
T~dt 



n 2 s = o. 



S{t = 0) = 5(0) = S x (0) + iS z (0) 



dS 



(26) 
0.(27) 



The solution of Eq. ( 26 ) with initial conditions ( 27 ) can 



be represented in the form 



S(t) = S(0)t 



cosh Kt 



sinh nt 
2tk 



(28) 



where k = -\/1/(4t 2 ) — Q 2 . It should be emphasized that 
the parameter k can take both real and imaginary values 
depending on system parameters. Moreover, in addition 
to describing spin polarization decay in infinite wires, Eq. 
( 28 1 also describes spin relaxation at long distances from 



boundaries of finite length (and semi-infinite) wires. 



Generally, we can distinguish two spin relaxation 
regimes. It follows from Eq. (28 1 that when 2t£1 < 1 



(Im(At) = 0), the spin relaxation is described by two ex- 
ponents with relaxation rates 1/(2t)±k. In the limit of a 
small spin precession angle per mean free path, rfi <C 1, 
we find 



S(t) = S(0)t 



rQ. z t 



(29) 



The above expression coincides with the spin relaxation 
time predicted by D'yakonov-Perel' theory for ID relax- 
ation. It is clearly seen that in this situation the relax- 
ation of spin polarization is characterized by the time 
constant (rfi 2 ) -1 , which is much longer than r. 

The second regime of spin relaxation is realized when 
2rf2 > 1. In this case k is purely imaginary meaning that 
the spin relaxation decay described by Eq. ( 28 1 consists 
of an exponential decay with a rate of l/(2r) modulated 
by oscillating functions. In the limiting case rJ7 > 1 we 
obtain 

S(t) = S(0)e- t/(2T) cosnt. (30) 
Both regimes of spin relaxation are presented in Fig. [I] 



Relaxation of inhomogeneous polarization in 
infinite wires 



As it is mentioned below Eq. (24), the Green's func- 



tion of Eq. (19), Sc(x,t), can be obtained performing a 
back transformation (from u to S) of the knowrP^l Green's 
function of Eq. ( 24 ) (the forward transformation is given 
by Eq. (23)). Following this procedure we find 



S G (x,t) = —$(vt-\ x \)e 
2v 



UjX 



.JL / V : v 2 t T 



2vt 



,(31) 



where $(£) is the Heaviside step function, and Io{t) is 
the modified Bessel function of zero order. We note that 



the Green's function (31) of Eq. (19) describes the evo- 



lution of a point excitation of spin polarization with the 
following initial conditions 



S G (x,t) = for t < 0, 



dSc 



dt 



5(x). (32) 



The Green's function (31 ) can be employed to determine 



the spin polarization at any moment of time for any given 
initial conditions using the relation 



S(x,t) 



d_ 1 

at + t 



oo 



-oo 
oo 
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FIG. 1: (Color online) Dynamics of spin relaxation of ho- 
mogeneous spin polarization in infinite wires, (a) and (b) 
correspond to two different regimes of spin relaxation of Eq. 
(28 1 as discussed in the text. The values of k are indicated 
on the plots. 



In order to better understand the meaning of Eq. (311, 
let us consider the spin dynamics of the following initial 
excitation of spin polarization in z direction: 



S(x,t = 0) = , 



dS(x,t) 



dt 



iS(x). (34) 



Accordingly to Eq. ( 33 1 , at any t > the spin polariza- 



tion in the system is given by 
S x {x,t) = 



<&(vt — \x\) _j_ I \Jv 2 t 2 — x 2 

smrjxe ^ I 



2v 



S z (x,t) = 



$(vt- \x\) 
2v 



cos rjxe 2t Iq 



2vt 
V 'v 2 t 2 — x 2 



2vt 



,(35) 



(36) 



This solution describes a propagation of the initial ex- 
citation of spin polarization in both directions from the 
excitation point with well defined fronts. Such a propaga- 
tion also involves the spin precession and spin relaxation. 
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FIG. 2: (Color online) Dynamics of S z (a) and S x (b) com- 
ponents of spin polarization of an initial Gaussian spin polar- 
ization profile with spin polarization pointing in z direction 
(S{x,t = 0) = iexp(-1002: 2 /^))- The spin polarization com- 
ponents are calculated using Eq. < |33[ ). This plot was obtained 
using the parameter value rjl — 0.1. 



Fig. [2]shows the dynamics of spin polarization in an in- 
finite wire found using Eq. ( 33 1 with different initial con- 
ditions. Specifically, we assume a Gaussian spin polariza- 



tion profile and zero time derivative of spin polarization 
at t — 0. This plot shows that the initial spin polarization 
profile splits into left- and right-moving packets whose 
amplitude in z direction decreases in time. At the same 
time, the amplitude of 5a, in packets increases because of 
the spin precession. It should be emphasized that such 
moving packets of spin polarization are not c aptured by 
the drift-diffusion schemes of spin transport! ^ 8 * 32 ! 33 !. In 
addition, we note that an area with a flat S z establishes 
between the moving packets. This region of flat distri- 
bution of S z can be considered as a precursor of diffusive 
dynamics. 



Moreover, the expression for the Green's function (31 ) 
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takes a simpler form at some specific points. For example, 
if we consider the position of the moving front, then it is 
easy to obtain 



S G {x = vt-Q,t) 



1 



(37) 



This expression basically means that the amplitude of 
the moving front exponentially decreases with the time 
constant 2t. Another interesting point is x — 0. We find 
that at x — the Green's function (31 ) is a real function 
of time: 



(38) 



Its asymptotic values at short and long times can be cal- 
culated. In particular, at short times Kt 



S G (0,t) 



I f I / 



It 2 /f 3 

o 1 



2v y 2t 16 t 2 
and at long times t 3> t 



a °M = s?Vs( 1+0 (i 



(39) 



(40) 



which corresponds to the well-known diffusive behavior. 



C. Relaxation of homogeneous polarization in 
finite length wires 
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Let us consider the problem of spin relaxation in wires 
of finite length L, < x < L, when at the initial moment 
of time the electron spins are homogeneously polarized 
along z-axis, 



5(x, 0) = iSo 



as 

~dt 



0. 



(41) 



We find the solution of Eq. ( 19 1 with the boundary con- 



ditions ( 20 1 and initial conditions ( 41 1 using the standard 
method of separation of variables. i\ straightforward ap- 
plication of this method leads to the following expression 
for the complex spin polarization 

■Sfo*) = ^(vL/ 2 ) c in(x-L/2) + 2„Le inx - t/(2T) 
S r)L/2 

^ 1 - (-l) ,l e-" fi / , sinh/c„A nnx ftn . 

x 22 r-T V. /--w ( cosh K "t + -^r^- ) cos — ,(42) 



2>TK n 



L 



n 2 n 2 v 2 1 L 2 



where n n = ^ 1/ [4rr" 

Previously, we investigated the spin relaxation in fi- 
nite length wires in the diffusive regime and found that 
an initially homogeneous spin polarization in z direction 
transforms into a persistent spin helix at long times^. 
The same qualitative result follows from Eq. (42 1: note 



that only the first term in the right-hand side of Eq. 



FIG. 3: (Color online) Dynamics of formation of persistent 
spin helix in the ballistic regime of spin transport. This plot 
was obtained using the parameter values r/l = 2 and L = 21. 



(42 1 survives at long times. However, Eq. (42) by itself 
is more complex compared to those reported in Ref. [S3] 
for the case of diffusive spin transport because it incor- 
porates both diffusive an d ballistic spin transport. It can 
be shown that Eq. ( [42] ) gives the same resuliP^l in the 
diffusive limit. Indeed, when £ = vt <C L and t r it is 
easy to demonstrate that 



lim e' t/(2r) { cosh K n t 

-s-CU-s-0 \ 



sinh n n t 



,(43) 



where D = £ 2 /t. Substituting Eq. Q into Eq. (| 42| ) 
we obtain the same expression for S(x, t) as in Ref. 1281 
The dynamics of formation of persistent spin helix in 
the ballistic regime of spin transport is depicted in Fig. 
[3] This fi gure is obtained plotting Eq. ( |42[ ) at different 
moments of time. One notices that the persistent spin 
helix is almost entirely formed within a short ~ 5r time 
interval. 

Let us consider the ballistic limit in more details. In 
this case assuming £ ^> L, 2k„t/z ~ 2irn£/L ^> 1 the 
second term in the big round brackets in Eq. ( 42 ) can be 
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omitted (since it involves a small multiplier 1/(2k„t)). Then, we obtain 



S(x,t) _ sm(r)L/2) in(x _ L/2) 



+ 00 



So 



V L/2 



L i V x-t/{2r) \- 1 - (-l)"e lr)L 



■nnix + vt) Tm(x 
cos h cos - 



vt) 



L 



(44) 



We can sum the series in Eq. ( 44 ) taking into account the 



fact that at the initial moment of time, t = 0, the right- 
hand side of Eq. (44) is equal to i. Moreover, since these 



Fourier series are even and 2L-periodic, we can finally 
obtain a very simple expression for the spin polarization 
in the ballistic limit: 

gOM) _ , sin(7?L/2) iTl(x - L /2) (, _ -t/(2r)\ 
So vL/2 V J 

where z = z + 2Ln, and n is an integer number selected 
in such a way that for any real z, —L < z < L. 



D. Relaxation of inhomogeneous polarization in 
finite length wires 

Next, we investigate the dynamics of relaxation of an 
arbitrary spin polarization profile in finite length wires. 
For this purpose, we obtain an expression for the Green's 
function G(x,t;£) that allows finding the spin polariza- 
tion profile at any moment of time for given initial con- 
ditions. The Green's function G(x, t; £) is calculated by 
mapping Eq. ( 23 ) into Eq. ( 24 ) taking into account the 



boundary conditions ( 25 1 by an even mapping of the ini- 



tial condition on an infinite wire. It is straightforward to 
show that G(x, t; £) can be written as 



G(x,t;0 



,iTi(x-£)-t/(2 T ) +°° 



2v 



Mvt-\x-t-2Ln\) 



( ^/vH 2 - (x - £ - 2Lnf 

^0 ^ 

\ 2vt 



<b(vt- \x + £-2Ln 



, T ( ^JvH 2 -{x + Z- 2Lnf 

> 1q o 

\ 2vt 



The spin polarization distribution along the wire at any 
moment of time can be calculated using the following 
equation: 



S{x,t) = 



d_ 
Ft 



+ / deG(aj,t;OS(f,0). (47) 
Jo 



As an example of spin dynamics calculation in finite 
length wires, let us consider the evolution of a Gaussian 
initial spin polarization profile. The time derivative of 
spin polarization at t — is set equal to zero (this is the 
second initial condition). Fig. [4] shows results of our cal- 
culations based on Eqs. po) and ( [47] ). Specifically, it is 
clearly seen that the initial spin polarization profile splits 
into two packets moving in opposite direction similarly 
to the result shown in Fig. [2] for the case of infinite wire. 
The amplitude of S z in these packets decreases in time. 
At t = 2t the left-moving packet reaches the wire bound- 
ary and is reflected back changing its direction of motion. 
We call such a reflected packet of spin polarization as the. 
boundary spin echo to distinguish it from the spin echo 
in nuclear magnetic resonance (NMR^pS. The boundary 
spin echo is indeed much closer reassembles the echo of 
sound waves than those in NMR. We also note that at 
long times any initial spin polarization profile (including 
one shown in Fig. [4]) shapes into the persistent spin helix 
characterized by a certain amplitude and phase. 



IV. SPIN RELAXATION IN 2D CHANNELS 

Let us consider a two-dimensional electron gas in the 
presence of linear in the wave vector Rashba and Dres- 
sclhaus spin-orbit couplings. The electron Hamiltonian 
is written as 



H = — + a {p y cr x - p x a y ) + P{p x a x 



PyCJy), (48) 



where a and /? are strengths of Rashba and Dresselhauss 
spin-orbit couplings correspondingly. In the specific case 
that we consider below, namely when a — /3, the Hamil- 
tonian (48) takes the form 



.(46) 



H = h a(a x - <r y )(p x +p y ). 

2m 



(49) 



Next, consider a rotation about z axis by an angle ip = 
7r/4. Then, x is transformed into x' = (x + y)/y/2 and 
y' = (y — x)/y2. The transformation of the spin part 
of the Hamiltonian ( 49 1 is performed using the standard 



unitary transformation 



U 



e -iy/ 2 



(50) 
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FIG. 4: (Color online) The boundary spin echo in a finite 
length wire. The evolution of an initially Gaussian spin 
polarization profile pointing in z direction (Six, t = 0) = 



iexp(-100a; 2 /^)) 

was found using Eqs. mm and &47«. It is 

clearly seen that the spin polarization packet initially mov- 
ing in —z direction moves back after the reflection from the 
boundary. This plot was obtained using the parameter values 
r)l = 0.1, L = 6L The curves are displaced for clarity. 



In the new basis the Hamiltonian (|49| takes the form 
H' = UHU ] = 



2m 



2ap x >cr ',, 



(51) 



where a' y , = a y corresponds to the projection of spin 
operator in the rotated basis on y' . Under the action 
of the Hamiltonian H', the spin of an electron precesses 
with an angular velocity tt' p — —fl'e y i cos 6*, where 9 is 
the angle between the direction of electron momentum 
and x' axis and f2' = Aap/h. 

Now, let us consider electron spin dynamics in a chan- 
nel of width L in x' direction and of infinite length in 
y' direction. Our goal is to show that the evolution of 
spin polarization (homogeneous in y' direction) in such a 



channel is similar to the evolution of spin polarization in 



finite length wires considered in Sec. Ill Assuming that 
the spin polarization in the channel depends only on x' 
and S y i — 0, we can combine Eqs. ^ and ^ as 



aS(fl) 
dt 



V cosf 



dx' 



-fi'cosfle^ x S(0) 



d0 

2^T 



8(6) 



(52) 



where the following shorthand notation is used: S(8) = 
S(6,x',t). Introducing a complex spin polarization 
S(9) = S x '(8)+iS z (9), we arrive to the following integral 
equation 



dS(0) 
dt 



+ v cos i 

1 
r 



} dS(9) 
dx 

S(6)- 



ift cos 9S(9) 



d9 
2^ 



S(6) 



complimented by the boundary condition 
(S(6)-S(7r-6))\ r = 0. 



(53) 



(54) 



We seek the solution of Eq. ([53]) in the form of Fourier 
series 



+oo 



S(6) = (Tq + > (<r„ cos n9 + a n sin n9) , 



n=l 



(55) 



where a n and a n are functions of x' and t. In particular, 
it is clear that a = d9 S (9) / (2n) is proportional to 
the spin polarization. 

Substituting Eq. ( [55] ) into Eq. ( [53] ) and using the 
completeness of the set of function {1, cos nd, sinn^} we 
find infinite series of interconnected equations. The two 
first equations in the series are: 



~dt 



dap 
dt 

d 



1 / d_ 

2 \ dx' 



V d^- m ' 



cto 



ifl'J o x = 0, 
- \ ■ — £ = o. 



(56) 
(57) 



For a weakly anisotropic electron momentum distribu- 
tion I (72 1 <C I Co |, therefore, we can truncate the series 



by setting (T 2 = in Eq. (57). Thus, Eq. (57) can be 
rewritten as 



da\ 
~dt 



— 
'dx 1 



in' ) (j + — = o 



(58) 



Excluding o\ from Eqs. (|56|), (|58|, we obtain the follow- 
pat io 

d 2 a 



ing equation for the total spin polarization Co 

1 / V d ' V 1 ' 

^dx 1 



dt 2 



1 da 
t dt 



n 1 



V2 



(Jo = (59) 



which is of the similar form with Eq.(19). This is the 
main result of Sec. |IV| proving similarity of spin relax- 
ation in finite length wires and 2D channels. 
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The boundary condition for spin dynamics in the chan- 
nel can be obtained in the following way. Substituting the 
approximate solution of Eq. ( 53 ) , 



S(9) = ctq + o\ cos 9 + (j\ sin 9, 



(60) 



into the boundary condition Eq. (54), we immediately 



notice that the function o\ on the channel boundary must 
turn to zero, namely 



o-i| r = 0. 



(61) 



Combining Eq. ( 61 ) with ( 58 1 , we find the corresponding 



boundary condition for the complex spin polarization gq, 

= 0, (62) 



f da ° 
I cV 



l (T 



which is similar to the boundary condition in finite length 
wires given by Eq. ( 20 1 . 



V. CONCLUSIONS 

In this paper, we studied spin relaxation in ID and 2D 
systems with spin-orbit interaction. For this purpose, we 
used the spin kinetic equation that takes into account 
both ballistic and diffusive spin transport regimes. This 
formalism was applied to several interesting problems of 
spin relaxation that were solved analytically. In partic- 
ular, we analyzed dynamics of homogeneous and inho- 
mogeneous spin polarizations in infinite and finite length 



ID wires. Moreover, it was explicitly demonstrated that 
the problem of spin relaxation in appropriately oriented 
2D channels with Rashba and Dresselhaus interactions 
of equal strength can be reduced to the problem of spin 
relaxation in ID wires. This result establishes a solid 
foundation for a recently suggested method of creation 
of persistent spin helical configurations in semiconduc- 
tors^. 

One of the main results of the paper is a prediction of 
the spin echo in the ballistic regime of spin transport. We 
found that in finite length wires an initially localized spin 
polarization profile reflects from a sample boundary re- 
turning to its initial position. To distinguish such a spin 
echo from that in NMR and emphasize its closer analogy 
with the spin echo of sound waves we named such an ef- 
fect as "the boundary spin echo" . Another interesting re- 
sult is a discovery of a transformation that maps the spin 
kinetic equation into the Klein-Gordon equation with an 
imaginary mass. In the relativistic quantum mechan- 
ics, this equation is used to describe tachyons, namely, 
hypothetical subatomic particles that moves faster than 
light. Therefore, we believe that certain predictions of 
relativistic quantum mechanics can be laboratory tested. 
Of course, much work is to be done in this direction. 
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